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Abstract 

The reaction rate probability integral is extended from Maxwell-Boltzmann approach to a 
more general approach by using the pathway model introduced by Mathai [Mathai A.M.:2005, 
A pathway to matrix-variate gamma and normal densities, Linear Algebra and Its Applica- 
tions, 396, 317-328]. The extended thermonuclear reaction rate is obtained in closed form 
via a Meijer's G-function and the so obtained G-function is represented as a solution of a 
homogeneous linear differential equation. A physical model for the hydrodynamical process 
in a fusion plasma compressed and laser-driven spherical shock wave is used for evaluating the 
fusion energy integral by integrating the extended thermonuclear reaction rate integral over 
the temperature. The result obtained is compared with the standard fusion yield obtained 
by Haubold and John in 1981. [Haubold, H.J. and John, R.W.:1981, Analytical representation 
of the thermonuclear reaction rate and fusion energy production in a spherical plasma shock 
wave, Plasma Physics, 23, 399-411]. An interpretation for the pathway parameter is also 
given. 

Keywords: Thermonuclear function, Mathai pathway model, Tsallis statistics. Guderley 
model, fusion energy, Meijer G-function. 

1 Introduction 

Energy production in the Sun and other stars is mainly due to thermonuclear re- 
actions taking place in them [9, 6, 2J. Hence the nuclear reactions in the laboratory 
situation for energy production has got more importance especially when the reactions 
are laser-induced fusion reactions [3j [T]. The fusion reactions are controlled by the 
thermonuclear reaction rates under specific theoretical assumptions and experimental 
inputs. The major facts which are required for the calculation of the reaction rates are 
(1) the quantum-mechanical cross section of the reaction and (2) the statistical mechan- 
ics velocity distribution of the reacting particles. A systematic study of reaction rate 
theory has been done by many authors [2H [19], [16] . In the literature it is assumed that 
the physical parameters in the solar interior (particularly density and temperature) lead 
to an equilibrium velocity distribution for fusion plasma particles. This distribution of 
particles is assumed to be Maxwell-Boltzmannian in almost all cases in stellar physics 
and cosmology [18] . A small deviation to the distribution of the reacting particles re- 
sults in a change in the reaction rates. It was observed that the distribution function 
in the solar interiors cannot be much different from Maxwell-Boltzmannian [U [7] . 
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The dynamics of an imploding spherical shock front in inertial confinement fusion 
was first investigated by Guderley using a self-similar form of the fluid equations [12]. 
This analysis was continued by others [2T| IT5] and has been reviewed extensively by 
Zel'dovich and Raizer [39J and Atzeni and Meyer- Ter-Vehn [I]. Brueckner and Jorna 
[3], Haubold and John [T7], and Rygg citerygg2006 calculated analytic approximations 
to the respective fusion energy yields. In such calculations, Maxwell-Boltzmannian ve- 
locity distributions are assumed for reacting particles but most likely this assumption 
is violated due to non-locality of the constituent particles and the convergence of the 
fuel at velocities comparable to the typical ion thermal velocity [32J. 



The spatio-temporal dynamics of solar activity has been investigated by studying 
the main solar oscillations and the time evolution of the basic periods [35]. Period 
length variations could improve knowledge of the relationship between observed solar 
activity due to variable phenomena generated by solar dynamo models or even in the 
gravitationally stabilized solar fusion reactor [361 EH EH [11] due to effects of a deep 
mixed shell on solar g-modes, p-modes, and neutrino flux. 



Tsallis, in 1988, suggested a generalization of Boltzmann-Gibbs statistical mechanics 
with a ^-exponential function as distribution function and q — > 1 recovers Boltzmann- 
Gibbs statistics [311 E] • An attempt was made to extend the theory of nuclear reaction 
rates from the Maxwell-Boltzmann approach to the Tsallis approach by many authors 
[331 EH ESI ESI 122] • In 2005, Mathai (231 HI] discovered a pathway model which covers 
the Tsallis distribution and many other distribution functions as special cases. Mathai 
[23] mainly deals with rectangular matrix-variate distributions and the scalar case is a 
particular one. The real scalar case of the pathway model has the following forms: 

fi(x) = ci|x| 7 [l-a(l-a)|x|' 5 ]^, a > 0,6 > 0, l-a(l-a)|x|* > 0,7 > 0,a < 1 (1.1) 

and fi(x) = when 1 — a(l — a)^! 5 < 0, where a is the pathway parameter. When 
a < 1 the model is a generalized type-1 beta model in the real case. When a > 1 we 
have from (11.11) 

f 2 (x) = c 2 \x\' y [l + a(a - l)\x\ s }~^ , - oo < x < oo (1.2) 

which is a generalized type-2 beta model for real x > 0. When a — y 1 the above two 
forms will reduce to 

h{x) = c 3 \x\^e- a ^\ (1.3) 
where the normalizing constants are given by 

c\ = -. t -i s t~* s , a < 1, (1.4) 

r(^)r(^ + i) 

^^5P?)' S M-T j>0,a>l (1.5) 



2 



and 



8{arf) & 



a -> 1 



(1.6) 



c 3 = 




where 3?(.) denotes the real part of (.). For different special cases of the above forms 
see [23, 28J. The pathway model can also be established from a generalized entropy of 
order a using the maximum entropy principle [231 [25] . Haubold and Kumar in 2008 
[18] extended the reaction rate theory to a general form covering the Tsallis reaction 
rates and established the closed form representation of the extended thermonuclear re- 
action rate in terms of generalized special functions, mainly in Meijer's G-function and 
//-function [22] ED3 ED 122] • 

In this paper we develop a complete theory for the representation of the nuclear 
fusion yield integral by using the closed form representation of the extended thermonu- 
clear reaction rate. Brueckner and Jorna [3] give a physical foundation of laser-induced 
nucleosynthesis. Instead of using difficult physical descriptions of laser-induced fusion, 
a simple physical model of Guderley [T2] has been used and fusion energy yield integral 
has been evaluated in closed form via a G-function. From the literature it is observed 
that Maxwell-Boltzmann distribution function is used in the model calculations con- 
cerning fusion reactions in strong plasma shock waves. But we use pathway model 
[23~t [28] which has a more general distribution function than Maxwell-Boltzmann distri- 
bution, for the calculations. As the pathway parameter a — >■ 1 the Maxwell-Boltzmann 
distribution is retrieved. 

The paper is organized as follow: In the next section we extend the thermonuclear 
reaction rate using the pathway density function and present the result in closed form 
in terms of a G-function. A homogeneous linear differential equation which satisfies the 
G-function is also obtained. In section 3 we illustrate the model formation by using the 
Guderley model [12] for the compressed fusion plasma and heated by a strong spherical 
shock wave used by Haubold and John [17] and Brueckner and Jorna )3] . Section 4 
gives the exact evaluation of the fusion yield integral by the extended reaction rate. 
A comparison of the new result obtained with the result already obtained by Haubold 
and John [T7] is done in section 5. A physical interpretation for the pathway parameter 
in the newly obtained result is also given. Concluding remarks are added in section 6. 

2 Evaluation of the extended thermonuclear reaction rate in 

closed form 

We consider the reaction between two interacting particles of certain type % and j, 
then the reaction rate is the product of the number densities rii and rij and the 
reaction probability (crv)ij between the particles which is given by 



The bracketed quantity (<Jv)ij is the probability per unit time that two particles i and 
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j confined to a unit volume will react with each other. It is the statistical average of 
the reaction cross section a(v) over the normalized distribution function n(v) of the 
relative velocity v of the center of mass system: 

POO 

(crv)ij — dv n(y)a{v)v. (2.2) 
Jo 

When we consider the cross-section a(E) for low-energy nuclear reactions far from 
any resonance, we obtain the expression for the cross-section of the low-energy non- 
resonant nuclear reactions as [HJ [9] 

a(E) = ^e-^\ (§«l) (2.3) 

where E is the relative kinetic energy between the particles, B the nuclear barrier height 
and e~ 2nv ^ is usually called the Gamow factor, gives the probability that the particles 
penetrate the potential wall. The Sommerfeld parameter 1](E) is given by 

V(E) = (^Y^ (2.4) 
v z ' TiE 2 

where Zi and Zj are the atomic numbers of the nuclei i and j, e is the quantum of 
electric charge, h is the Planck's quantum of action. The cross-section factor S(E) 
entering the cross-section is a slowly varying function of energy over a limited energy 
range and can be characterized depending on the nuclear reaction. Usually when 5" is a 
constant one takes it as So- S(E) contains the constants intrinsic to the nuclear reactions 
under consideration. Its energy dependence may be represented by a Maclaurin series 
expansion up to the second order in the kinetic energy E as following [T6| ITT] : 

S(B) „ S (0) + ^ + i^ = t^. P.5) 

Due to the strong dependence of the cross-section a(E) on the relative kinetic energy 
E at the temperature of thermonuclear burning, the essential contribution to the ther- 
monuclear reaction comes from the collisions between particles with energies more than 
the average thermal energy (^r). Hence the reaction probability depends on the re- 
acting particles with energies greater than the mean thermal energy. 

Usually, the thermonuclear fusion plasma is assumed to be in thermo dynamical 
equilibrium with regard to their velocities (not with regard to their mass abundance). 
From the tabulations of the nuclear reaction rates by Fowler et al. [10] , we can infer 
that the distribution of the relative velocities of the reacting particles always remains 
Maxwell-Boltzmannian. In case of non-degenerate and non-relativistic gas the distri- 
bution function of the relative velocities of the nuclei is Maxwell-Boltzmannian. 

3 

fMBD(E) = 2tt f-jLj 2 e-^v^EdE. (2.6) 
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By using the Maxwell-Boltzmann distribution function (12. 6ft . the non- resonant low- 
energy cross-section (12.31) and the expression for the cross-section factor given in (I2.5p . 
the reaction probability in (I2.2p becomes 



(av) 



i.) 



TTfl 



1 
kf 



E 

i/=0 



E"e K J ^^dE. (2.7) 



If we assume a deviation from the Maxwell-Boltzmann distribution function or a more 
general distribution function, then we can consider a distribution function fMPD(E) 
obtained by the pathway model of Mathai (23J, [25] given by 



fuPD(E) = 2ix- 



E 



irkT) 



1 + (a 



E ' 
'kf 



i 

' a-l 



dE. 



{21 



As a —7- 1 in (12.81) we obtain the Maxwell-Boltzmann distribution function given in 
(12. 6p . Substituting Mathai's pathway distribution (12. 8p instead of Maxwell-Boltzmann 
distribution (12. 6 p in ( 12. 2p we get, 



7T/X ) 1 fcT ) S 



5— S'M(O) 



y=0 



E" 



! + («»- l4 



1 

' a-l 



e V2; dE. 



(2.9) 

For the convenience in evaluating the above integral we give a substitution y = and 



x 



My = 



2 ZiZje^ 



7T/i 



I 2 

2 



. Then the equation (12. 9p becomes 



E 

i/=0 



1 



+ (a - l)y]-^ie- a,tf h dy. (2.10) 



The reaction probability (crv)ij is the finite sum of integrals of the form 



ha 



y"[l + (a-l] 



! dy. 



(2.11) 



Using the Mellin convolution property the above integral can be evaluated and expressed 
in closed form by using generalized special functions. The function can be expressed 
in terms of a Meijer's G-function of one variable. G-function G'^^ l (z\^'^'"'^) was 
introduced by Meijer in 1936 as a generalization of the hypergeometric function, see 
[22| |30~1 129] for details. For the detailed evaluation of ha in (12. lip see [18]. Thus we 
obtain 



ha 



(tt)-I 



-Gli 



a 



l)x 2 



1 0,|,H-1 



(« - i)^ l r (A) 

According to this the extended reaction probability in (12. 9 j) becomes 

/J_\ _I/+ ^ 5-^(0) 1 

i/=0 v 7 



(2.12) 



{crv)i 



xGli 



2kT) V h 



lo,i^+i 



(2.13) 
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The Meijer's G-function G^g ^z\ j satisfies the homogeneous linear differential 

equation of order three [30j 129] • 

3 



i 2 1 - T +u 







- : I z— — ai + 1 



-z[ z— + 



dz a — 1 



-n 

i=i 
1 - v 



z Tz - b > 



1=0,62=5,63=^+1 



_d_ / _d _ 1 

Z dz \ dz 2 



x z 



dz 
dz 3 



l + u 

'3 



/-<3,1 I \~ a-l ' 

G l,3 I « I o.I.H-1 



2 d 2 /z/ \ d 

^ I 2 -rrr + z) z— 

2 / dz 2 V2 / dz 



+ 1 + 1/- 



1 



a: 



/^3,l 1 1- a-l 1 
^1,3 I ^lol^+l 



(2.14) 



Substituting 2; 



x 



d*_ 

dx 3 



2u 



dx 2 



(a-l)x 2 



+ (a 



in (I2.14p we obtain 



l)x— -2[(1 + v){ a - 1)-1] 



3,1 
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l)x 2 ,2— 
I o,i 



-+1/ 



(2.15) 

3 Fusion yield integral in the case of Shock-compressed and 
heated plasma 

To produce a useful amount of fusion energy the plasma pellet is to be highly and 
efficiently compressed. The center of the compressed pellet must be brought to the 
ignition temperature and density, but the rest of the compressed pellet must be left 
as cold as possible. Very high compression is needed to maximize the reaction rate 
after the pellet is ignited [3j. A laser can be used as a hydrodynamic driver of pellet 
compression. The uniform pressure applied to the pellet surface by the laser energy 
will produce a single converging shock [3 |3]. We take the Guderley model (12] of 
the dynamics of strong spherical shocks near the center of the sphere provided the 
perturbation of the hydrodynamic process by the fusion reaction is ignored. As the 
converging shock approaches the center of collapse, it decouples from its generating 
boundary. In a strong spherically convergent shock the shock position r can be written 
as a function of time to as 

r = a-to) n , (3.1) 

where we measure the time at which the shock reaches the center of convergence. The 
parameter £ measures the shock strength, for the number n, < n < 1 [12]. In spherical 
geometry and for an ideal gas with the ratio of specific heats 7 = |, the exponent n 
has to take the value n 
equation (13.1 1) as 



0.688377. The shock velocity can be determined from the 



u 



dr 

dF 



n£(-it 



\n— 1 



t- 1 



(3.2) 
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where the time to is negative and it increases up to the center of convergence at to = 0. 

The scaled density, scaled temperature and scaled pressure at a radius r are functions 
of the reduced time — where — tn is the time after which the shock front has reached 

to u 

the center of convergence. 

P f,(-f), (3.3) 



Po V h 

— 2 =f P (~l) (3- J) 



mp muo 2 fpy~h) mu 2 ( t 



- 2, 2 /pH) 2 <X5) 

where m is the mean mass of ions, electron and ion temperature are assumed as equal, 
6 is written for kT. Also fg = jr. For evaluating the functions in 03.3p .( l3Tl]) and (13. 5p 
for density, pressure and temperature the initial conditions can be taken as the jump 
conditions. When the first converging shock is passed a discontinuity is produced at 
the center. The computations of the functions f p ,f p and fg can be seen in Goldman 
[14J. The fusion energy can be evaluated by integrating the thermonuclear reaction rate 
depending on position and time via density and temperature. 



r^max r*max 

E 'fusion = 47r Eij r 2 dr rij(r,t)dt (3.6) 

where denotes the energy released by a fusion reaction between the nuclei i and j. 
Using the reaction rate given in ( 12. ip we have 

Tij = ( l -\h) N i N M v )ii 
1 iV 2 

= (l + 5 ij )- 1 — {av) ir (3.7) 

Here N denotes the total ion density of the fusion plasma, 5{j is the Kronecker delta. 
The factor | is due to the fact that each reactant density is half the total ion number 
density. Now take N = where iV is the ion density and p is the local density. The 
dimensionless variable r = — ^ which is positive since t < 0. We have 



E fusion = tt(1 + <),,) ] \rE,, J r 2 (J J J f P \r){av = F{9{r, r)))ydrdr 

(3.8) 

where A^ = f ,r max = -W. 

From Goldman [13] we can see that the passage of the first converging shock gives 
a density increase followed by an adiabatic compression to a density ratio. Then the 



i_ 

max / rr> \ n /" ' max 
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reflected shock at the center on returning gives a further shock compression to a maxi- 
mum density [3]. It should be noted that temperature is a rapidly varying function of 
r but slowly varying function of r in the reflected shock. So in ( 13. 8 p we can write 

(<ru = F(9(r, r)) y = (av = F(9(r, n))}^- (3.9) 

where T\ is the scaled time of the reflected shock. So in the integral with respect to r 
we can take Ti = — j- instead of —1 = For smaller values of r and near the time 

to to 

of maximum compression the fusion yield Ef US i on shows a rapid increase. So we can set 
the upper limits of integration in ( 13. 8 p as oo. Now (13. 8p can be modified as 

E fus%on = tt(1 + Si^NjEa Qj " J r 2+ ^ J f p 2 (r)(av = F(6(r, r))) ii drdr 

(3.10) 

In front of the secondary shock the relation connecting the density, temperature and 
pressure at a radius r is given by 

f = f P (r), (3.H) 

Ps 

^ = f P (r), (3.12) 

Ps 

mp mp s f p (r) mp s 
2p 2p s f p {r) 2p s 

The equations are from Haubold and John [17]. By the jump conditions across the 
shock [El Q3] 

p s = Po 1 ^,l^±l (3.14) 
7 - 1 

7 + 1 

and using equation (13.21) the temperature following the secondary shock is a function 
of r given by 

a, x ™u Q 2 7 - 1 
7 + 1 7 + 1 

„ 2 2(n-l) 

run qnr « 7 — 1 
7+1 7+1 
From (13.161) the radical variable r is obtained as 



fe(n). (3.16) 



mn 2 £n 7 — 1 ^ 
M T i, 



2(n-l) 

n 

[6>(r, Ti)] 



7 + 1 7 + 1' 

A-^T) [0( r>Tl )]2js=ij ) (3.17) 
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where 



A 



mn 2 ^» 7 — 1 



fe(n) 



(3.18) 



7 + 1 7 + 1' 

Changing the variable of integration from the radical variable r to the thermal energy 
variable 9, the equation (I3.10p becomes 



E 



fusion 



where 



3n + l 

02(n-l) 



'0 

"OO 



2 1 -n) 



7r(l + 5 lJ )- 1 iV % 

r 00 +3 

tt(1 + S^N^E^h^ / = F(9))ijd9, 

Jo 



f p 2 {T)(av = F{9)) ij drd9 



h 



f P 2 (r)dr 



n 



2(1 -n) 



3n + l 
-\2(l-n) 



(3.19) 



(3.20) 



(3.21) 



Now we evaluate the fusion energy released by the plasma substituting the extended 
thermonuclear reaction probability over temperature. 

4 Evaluation of the fusion energy integral by the extended 
thermonuclear reaction probability in closed form 

With the help of the closed form evaluation representation of the reaction probabil- 
ity obtained in section 3 we find the closed form of the fusion energy yield integral 
(13.191) in the previous section. The calculations are done by taking the cross section 
as a non-resonant case in the energy region. It is found in literature that in the lab- 
oratory experiments, for example reactions like D(d,n) 3 He and D(d,p) 3 H, we take 
non-resonant case. In equation (12. 5p of section 2 if = for v = 1 and v = 2. Then 
we obtain the fusion yield integral after inserting the reaction probability (I2.9P in ( 13 . 1 9[) 
by writing 9 for kT. 



E 



fusion 



■ 1 JV 2 £7«/iC n S'(0) 



i)r(sb) 



la 



2 

— 



xGl 



3.1 

,3 



[a 



d9. 



By using the substitution 9 = -, we get 



E 



fusion 



i + ( y- 1 iV ^0^(0) 



("-i)r(^) 



(a - 1)tt 2 



2)\ h 



2\ 2 



u\ 



du. 



2n 

u 1 ~ n 
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By taking 



we get, 



E 



fusion 



h 



\fij (a-l)T(^) J 



x 



(4.1) 



Ml- 



3> 3 



4 k|,o 



dtt. 



(4.2) 



In view of the Mellin transform of a G-function, equation (14. 2 j) becomes 



/nsion ;: (1 + ^ij) ^0 EijIiC n S(0) , , . , 

>/ («- 1 ) r ( — ) 



(a — l)x 



v.2 



x rii + ?±iWi + ?±±W?±±V 1 



1 — n J \ 2 1 — n 



where a >l,^ T -l-S±i>0. 

' a— 1 1—n 

By Legendre's Duplication formula [22j 129] . 



r(2z) = 7rh 2z - 1 r(z)r(^ + ^ 



1 — n I \a — 1 



1 - 



71 + 1 

1 — n 



we get, 

E fusion 



(i + ■y-W-Eis/iG.sm (^ i 2 2-t 



1 



where 



1 — n 



1 — n 



(«-i)r(ir) 

2 \ 1 



(a — l)x 2 



(4.3) 



(4.4) 



n+l 



a — 1 1 — n J 1 a — 1 1 — n 



>0 (4.5) 



/pV)dr 



71 , 3n+l ( 1 

\2(l-n) 



2(1 -71 

2(l + n) 
71 1 ~ n . 4 

2(1 -nf 



m 7 — 1 



fe(ri) 



_, l + 3n 

2(l-n) 



_7 + 1 7 + r 

For a physically realizable solution we take n = 0.68837. Then equation (14.51) becomes 

1 



E 



fusion 



(l + 6 ij )- 1 N ' i E ij I 1 C n S(Q) 



X 



(«-i)r(^r) 

r(6.418)r(10.836)r - 6.418^ 



[/ t\ -21 -5-418 

(4.6) 
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Now we compare the newly obtained result with the existing results in the next section. 

5 Comparison of the results and interpretation of the pathway 
parameter 

The fusion yield integral in (14. 6p is obtained by using the extended thermonuclear 
reaction rate (I2.9p in (I3.19p . If we take the limit of Ef usion as a — > 1 we obtain the 
fusion yield integral in the Maxwell-Boltzmann case which is given in Haubold and 

John rm. 



EfLon = (1 + S^No'E^hCnSiO) (jy^j ' [x^ 5 ' 418 r(6.418)r(10.836) , (5.1) 

by using the asymptotic expansion of gamma function [8] 

T(z + a) ~ (2n)h z+a 'h~ z ,z ->• oo, \arg(z + a)\ < tt - e, e > (5.2) 
where the symbol ~ means asymptotically equivalent to. 

The interpretation of the pathway parameter used in extending the results can be 
done with the help of the 8 th moment. By using the type-2 beta form of the pathway 
model given in (11. 2p and the normalizing constant in (II. 5p . the (^-moment in the case 
of a > 1 is 

8[a(a - lM^T (+r) f°° ^ x i 

, x > 0, a > 1. (5.3) 



r(¥)r(^-^) 

7+1 

a[{28 + 7 + 1) - a( 7 + 1 + 8)} 



where £{.) denotes the expected value of (.). The value of a in terms of the 8 th moment 
is given by 

1 

a = 1 + 



(7 + 1 + 5) 



8- 7+1 



, x > 0, a > 1. (5.4) 



a£(x 5 ) 

For 7>0,5>0,a>0, the pathway parameter a increases and finally when the mean 
value or expected value of x 5 decreases a goes to — oo. While deriving the pathway 
model by optimizing Mathai's entropy one of the constraints was that £(x 5 ) is given. 
The value of a is determined depending upon the preselected value of £{x s ). When 
£(x s ) is set at high value then a approaches 1 + ^ + 5 l+S ) and if £(x 5 ) is a lower number, 

near to zero, then a will approach — oo. When a — > 1 we have £(x 5 ) = The values 
of £ (x s ) and a in the case of a < 1 is similar to (15 .3p and (15 .4p . Therefore 

< 1 c I 8\ < 7 + 1 

> ar]8 

In (I2.10p we have taken y = Here if T increases y will decrease and vice versa. 
If the expected energy £(^) < f then a < 1, if £(^) = | then a = 1 and when 
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> | then a > 1. The value of a in ( 12.81) is obtained by putting 7 = |, 8 — 1 and 
a = 1 in (15 .4p and is given by 

7 3 7 

a = t-ett =>■ a < - or a > —00. 

5 5£ (§) " 5 

6 Conclusion 

We extended the thermonuclear reaction rate from the Maxwell-Boltzmann approach 
to the Tsallis approach and obtained an analytical closed-form representation of the 
extended reaction rate probability integral. Instead of using the Maxwell-Boltzmann 
density function in equation (12. 6p we used pathway the density function (12. 8p and the 
extended reaction rate has been represented in terms of Meijer's G-function. We also 
obtained a linear homogeneous differential equation (I2.15P of order 3 satisfied by the 
G-function of the extended reaction rate. Then by using the Guderley model for the 
hydro dynamical compression and heating of the fusion plasma by a single laser-driven 
spherical shock, a general formula for the fusion yield integral is obtained by using the 
closed form representation of the extended thermonuclear reaction rate. The results 
obtained are compared with the result by Haubold and John [17]. An interpretation of 
the pathway parameter is also given. 
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